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We investigate the kinetics of normal phase nucleation and 
flux line condensation in the type-II superconductors by nu- 
merical study of the time-dependent Ginzburg-Landau equa- 
tion. We have shown that under the sufficient transport cur- 
rent the normal phase nucleates in the superconducting strips 
in a form of the macroscopic droplets having the multiple 
topological charge. We discuss the stability and the dynam- 
ics of the droplets. We found that pinning suppresses the 
droplet formation. 

PACS: 74.60.Ge,68.10.-m,05.60.+w 



The study of magnetic flux penetration in type II su- 
perconductors has attracted a wide interest both in view 
of important technological questions and as an proto- 
type of a general class of problems of nonlinear dynam- 
ics. Observations showed that the flux dynamics exhibits 
features that are similar to the viscous-fingering growth 
phenomenon in liquid-solid systems &-^|. In particu- 
lar recent experiments revealed dendrite flux penetration 
and the fingering of the remagnetization front. |^-|^] . The 
formation of vortex structure is traditionally viewed as 
the sequential penetration of vortices through the Bean- 
Livingstone surface barrier Q . It was found recently that 
flux penetration may also occur via dynamic instabilities 
of order parameter caused by the applied current and/or 
magnetic field. Numerical simulations revealed the inva- 
sion of the extended macroscopic normal areas (droplets) 
carrying flux into the superconducting sample [j^] . 

While the formation of normal areas looks natural for 
type I superconductors with the positive surface energy of 
the normal-superconductor (NS) interface, it seems sur- 
prising at the first sight that such an interface having 
in the static case the negative surface energy persists in 
type II superconductors. We see the explanation of this 
phenomenon in the fact that the transport current or al- 
ternating magnetic field, drives the superconductor into 
a strongly non-equilibrium state where the moving in- 
terface becomes stable. The idea that the free energy 
considerations do not apply to nonstationary processes 
in superconductors was put forward by Anderson et al 
Q in the context of phase-slips phenomenon. 

In this letter we report on our investigation of the ki- 
netics of normal phase nucleation and flux lines conden- 
sation in the type-II superconductors. We present the 



results of a numerical study of the dynamics of the flux 
penetration into strips with the transverse dimensions 
less than the effective penetration length X e ff = X 2 /h, 
where h is the thickness of the strip and A is the Lon- 
don penetration depth. We propose that the existence of 
the macroscopic normal regions is the direct consequence 
of their motion under the transport current. A current 
cannot penetrate the immobile compact normal zone im- 
mersed into a superconductor Q , therefore the NS inter- 
face moves towards normal phase with the velocity going 
to infinity |ic|ji~l)| , and the normal droplet disappears. 
At the same time, the penetration of the current into a 
normal phase makes it stable with respect to small fluc- 
tuations i.e. the transport current drives the system 
into a bistable state. Since the expulsion of the current 
from the normal regions requires a finite time, the cur- 
rent penetrates the moving normal droplet. The normal 
state develops and invades into the superconducting re- 
gion provided the current j flowing through the interface 
exceeds the stall current j* ]l^Jl^| . Thus the sufficient 
transport current stabilizes moving nuclei of the normal 
state in type II superconductors. 

The process of flux penetration occurs via the depres- 
sion of the order parameter on the macroscopic scale 
and can be viewed as the nucleation of the extended 
droplets of the normal phase in the superconducting sam- 
ple. An adequate description of such process involving 
the fast variations of the order parameter on the relevant 
spatial scale is given by the time-dependent Ginzburg- 
Landau equation (TDGLE) completed by the appropri- 
ate Maxwell equations: 
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where VP is the (complex) order parameter, ip — arg^ , 
A and \x are vector and scalar potentials, and j is the 
current density. The value of the dimcnsionless mate- 
rial parameter u is obtained from the microscopic the- 
ory The unit of length is the coherence length £, 
unit of time is to = £ 2 /Du, D = vfI/3 is the diffu- 
sion constant, I is a mean free path, vp is a Fermi ve- 
locity, the field is measured in units of the upper criti- 
cal field H c2 = <&o/27r£ 2 , $0 is the flux quantum. The 
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unit of current is jo = crTi/2eto, where a is the nor- 
mal conductivity. In these units the depairing current 
j p = 2/3-\/3 ps 0.3875. The condition X eff > 1 enables 
us to neglect the magnetic field created by currents 
and therefore drop the equation (^). We choose the ori- 
gin of the coordinate frame at the mid-point of the strip 
with the x axis lengthwise and the y axis in the lateral 
direction, so that the edges are located at (x, —d/2) and 
(x,d/2). The normal to the strip magnetic field B is 
associated with the vector potential A = (By, 0,0). 

We performed numerical simulations of TDGLE. We 
took a homogeneous superconducting state initial condi- 
tions (^ = 1, i.e. the state without magnetic field) per- 
turbed by a small amplitude noise. We used the no-flux 
boundary conditions, dy*£ = (i.e. the boundary with 
vacuum) in the transverse direction and the NS bound- 
ary conditions in the longitudinal direction (^(x, y) — > 
for x — > 0, L, where L is the strip length). We apply the 
split-step method described in fffijfy , the number of the 
grid points was 256x256 and the time-step was 0.05 — 0.1. 
Results of the simulations are shown in Fig. 1 . The sim- 
ulations were performed for j — 0.25, B = 0.0175 where 
as it had been shown in 0, the pure superconducting 
state is unstable with respect to vortex nucleation (note 
that our equations do not contain fluctuations). The in- 
tegration domain was 120 x 60. 

On the Fig. 1 the large dark droplets (for t — 50, 80) 
represent the normal phase emerging at the one side of 
the strip and traversing toward the opposite edge. The 
droplets are the long-living objects and as well as vor- 
tices play a crucial role in the dissipative processes. In 
our simulations the topological charge of these droplets 
would become as big as 5 -7 and even more. The droplets 
possess long tails (due to a finite relaxation time of the 
order parameter at the superconducting areas swept by 
the droplet). Our simulations show that new vortices 
appear at the edge just at the tail and then get sucked 
in the droplet. This can be easily understood since the 
formation of the new vortices is favored in the regions 
with suppressed order parameter. The normal phase ar- 
eas can evolve in two different ways. First, the normal 
droplet emerges at the edge, passes through the sample 
and vanishes at the opposite edge of the strip. In the 
second scenario that occurs under elevated currents the 
droplet traverses a strip leaving a channel (wake) of the 
normal phase behind. 

This scenario is shown on Fig. 1, t = 260. Than this 
channel traversing the sample breaks into the sequence of 
vortices (vortex street), which then propagate across the 
strip and annihilate at the edge. The nucleation and the 
propagation of the droplets and the vortices gives rise to 
non-periodic voltage oscillations along the strip. 

The droplets posses the topological charge n propor- 
tional to the gain in the superconducting phase along the 
loop enclosing the normal area. A relationship between 
the characteristic size R of the nucleus and n is deter- 



mined then from the condition that the supercurrent en- 
circling the nucleus (~ n/R) becomes equal to the j p giv- 
ing R ~ n/ j p . The size of the droplets in the strip can be 
estimated from the condition that total transport current 
at the distance R from the edge j(R) ~ j — B(R — d/2) 
is equal to the j p . It gives R = d/2 + (j — j p )/B. For 
the chosen parameters we obtain R ps 20 and n ps 6 -j- 7, 
which is in qualitative agreement with the results of sim- 
ulations. We expect that the above consideration holds 
also for the large, d 3> A samples where A takes the role 
of the characteristic length. We observed that droplets 
move much faster than single vortices. Simple analysis 
shows that the Magnus force exerted on the droplet grows 
linearly with n whereas the mobility saturates for large 
n, resulting in the velocity growth. 

Shown on Fig. 1 is a sequence of snapshots demon- 
strating a remagnetization process (we reversed the di- 
rection of the magnetic field at t = 200). At the first stage 
of remagnetization large normal phase areas develop at 
the edge of the strip. These areas swallow vortices corre- 
sponding to the previous direction of the magnetic field. 
Then the normal areas assume more complicated form 
and then break up into smaller droplets. For the zero 
applied current the Abrikosov vortex lattice is formed in 
the external field. In contrast to the case with nonzero 
applied current, vortices penetrate from both edges of the 
strip. When the direction of the field is reversed, large 
normal areas develop at both edges and swallow vortices 
corresponding to the initial direction of the field. After a 
while the new Abrikosov lattice forms with vortices along 
the reversed direction of the field. 

To include the Hall effect in our simulations we in- 
troduce the complex material parameter u — 5.79 + i. 
The imaginary correction to u describes the effect of the 
transverse Hall force on the vortices drift JItJ . This gives 
rise to the Hall voltage, moreover, we observe the turn 
of the droplet tail. We suggest that the rotation of the 
droplet 's tail in the experimental work || is caused by a 
significant Hall contribution. 

To summarize, we have found the long living droplets 
of normal phase inside a superconducting phase, and ob- 
served that they may posses the topological charge than 
can significantly exceed unity. Note that the droplets 
must be distinguished from the Abrikosov vortices with 
multiple charge. The linear stability analysis shows that 
such vortices are unstable with respect to the splitting 
into single charged vortices. The characteristic time 
of the splitting is of about 10—15 dimensionless units 
and, therefore, cannot explain the existence of long- living 
droplets. Note that these droplets may be viewed as the 
result of the "fusion" of the separate vortices. 

The qualitative arguments describing the droplet dy- 
namics can be put on the more rigorous basis for the 
droplets with the size well exceeding the coherence length 
£. In this case the boundary of the droplet can be con- 
sidered locally as slightly curved NS interface. Inside the 
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droplet the order parameter ^ vanishes and the field is 
described entirely by Laplace equation 

Afi = (5) 

The Eq. || has to be completed by the boundary condi- 
tions on the interface, deduced from the continuity equa- 
tion Vj = 0. It gives the relation between the com- 
ponents of currents normal to the interface jjfi = jn , 
where superscripts s,n denotes currents in normal and 
superconducting regions respectively. Using Eq. (||) 
we arrive at the first boundary condition — V n /z' n ) = 
|^| 2 (Vn<^— ■ A n ) — V„;j' s ' (here V„ means normal projec- 
tion of the gradient). The order parameter in supercon- 
ducting region near the slightly curved interface is given 
in "adiabatic approximation" by \^f\ 2 = 1 — (Vtpo — A) 2 . 

The phase ip of the superconducting order parameter 
in the leading order is described by Laplace equation 

= (6) 

together with the equation for normal velocity of the in- 
terface. They latter can be derived from Eq. (EI) for the 
slightly curved interface. The small curvature \ renor- 
malizes the normal velocity c n of the interface according 
Gibbs-Thomson condition c n = cq — \, where cq is the 
velocity of flat interface. 

For the flat NS interface the velocity Co(j) is a func- 
tion of the transport current. The one-dimensional sit- 
uation had been considered in [fl4|| , where the existence 
of the "stall" current j* at which the interface veloc- 
ity becomes equal to zero had been established. For 
u = 5.79 the stall current was found to be j* = 0.335, 
and co(j — ► j*) rs a(j* — j), where a — 0.6 is the numer- 
ical factor. In two dimensional situation the topological 
charge of the droplet induces the circular current j T tan- 
gential to the interface which modifies its the velocity. To 
account for the effect of the tangential current we take 
the order parameter close to the nearly flat interface in 
a form (the interface is parallel to y-axis, and we use 
a frame moving together with the interface along x-axis 
with the velocity c) — F(x — ct) exp[ik y y + <j)(x — ct)], 
where k x = ]ha x ^- 00 <f> x . and j T = (1 - k 2 - k 2 )k y , 
j„ = (1 — k 2 ~ ky)k x . A simple scaling analysis shows 
that the current renormalizes the interface velocity c as 

c(jn,3r)^C (j)^l-k 2 y , (7) 



where j — j n /(y 1 — k 2 ) 3 . If the curvature of the in- 
terface is small (i.e x — 1/R 1), the interface itself 
is defined by the additional condition that at the (flat) 
interface fi = (j,q = k x c(j n , j T ). After that the problem is 
completely defined. 

In the superconducting phase we have the Eq. (^) 
completed by the boundary conditions for tp on the strip 
edges. Thus the problem under study is a generalization 



of well-know problem of the Laplacian growth (see, e.g. 
Jl8|| ) . A new feature is that the function <p is a mul- 
tivalued one and has branch cuts. This multivalueness 
means that the obtained equations contain implicitly vor- 
tex solutions: vortices can appear and/or vanish via the 
formation of singularity at the interface. The detailed 
consideration of these equations we leave for the future, 
for now we would like to mention that the linear stabil- 
ity analysis shows that the flat interface with the cur- 
rent flowing through is stable with respect to small per- 
turbations. The above discussion and the results of our 
simulations make us to conclude that the passage of the 
current suppresses the NS interfaces instabilities in thin 
superconducting films. 

To study the effects of pinning we carried out sim- 
ulations of TDGLE with randomly distributed pinning 
centers. In the presence of the weak pinning the newly 
formed droplets assume the " fractal" configuration since 
the normal phase tries to settle at the pinning sites where 
the order parameter is already suppressed (see Fig. 2). 
The moving droplets percolates along the easy paths con- 
necting the pinning sites, but the pinning centers impede 
the interface motion. As a result the current that pene- 
trates the normal area gets smaller and cannot support 
the existence of the droplet any more and droplets break 
up. For stronger pinning the droplets do not form at all, 
and single vortices penetrate the strip via the jumping 
resembling the vortex motion through the array of linear 
defects g§. 

Finally we discuss briefly the time scale of the observed 
effects. The characteristic time in dirty superconductors 
is t Q ~ h/T c (l - T/T c ) « 10" 14 -7- 10~ n sec depending on 
the temperature interval. This means that the considered 
phenomena develop on the nanosecond scale. However 
the process of the flux penetration can be considerably 
slowed down by pinning. In this case the characteristic 
time (for "dendrite" formations for example) should in- 
clude macroscopic characteristics such as the size of the 
sample and the average pinning strength po| and can 
grow up considerably. 

In conclusion, we have shown that under the suffi- 
cient transport current the normal phase nucleates in 
the superconducting strips in a form of the macroscopic 
droplets which tear off at the edges and further propa- 
gate across the sample. These droplets possess the mul- 
tiple topological charge related to the magnetic flux they 
carry. Pinning suppresses the droplet formation convert- 
ing normal area into the multi-connected fractal forma- 
tions which then split into the separate vortices. We 
believe that the observed phenomena are not specific to 
the thin strips, and that the same mechanism governs the 
normal phase formation in large samples as well. 
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FIG. 1. Dynamics of the normal phase. The current is ap- 
plied along :r-axis and the magnetic field is perpendicular to 
the strip. Gray-coded images show ^(a;, y) | ( l^t = is shown 
in black and \^\ = 1 is shown in white). The field is reversed 
at t = 200. 



FIG. 2. Normal phase penetration at t = 40, j = 0.25, 
B = 0.018 in the presence of 180 randomly distributed pin- 
ning centers, other parameters are the same as in Fig. 1. 
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